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Abstract 



We have found in numerical simulations that the chaoticity of the classi- 
cal hamiltonian lattice SU(2) gauge field system is reduced in the presence 
I of static charges at the same total energy. The transition to non-chaotic 

' behavior is rather sudden at a critical charge strength. The equipartition 

. of chromoelectric and chromomagnetic energy takes place on a time scale 

essentially faster than the leading Lyapunov exponent. 

in 

&: 1 Introduction 



Knowledge of the mechanisms responsible for local equipartition of energy 
^ ' carried by non-abelian gauge fields and in particular knowledge of the speed 

^ by which the allowed phase space of possible field configurations is filled 

are important for the understanding of processes which lead to equilibrium 
in the very early universe and in relativistic heavy ion collisions. Prime 
examples for such processes are baryogenesis during the electroweak phase 
transition, the creation of primordial fiuctuations in the density of galax- 
ies, and the thermalization of excited quark matter in relativistic heavy 
ion collisions. Deterministic chaos of extended hamiltonian systems, like 
lattice models of non-abelian gauge theory are, is a universal mechanism 
of ergodization due to soft classical fields. Its main effect, the heating up 
of the system, is in addition to hard parton scattering. 

In the recent past numerical studies of the real-time evolution of classi- 
cal gauge field systems demonstrated chaotic behavior of randomly chosen 
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configurations with chromomagnetic energy The maximal Lyapunov 
exponent (measured in lattice units) has been found to scale linearly with 
the energy of the system (also in appropriate units) making thus an extrap- 
olation to the infinite dimensional field theory in the classical continuum 
limit possible. The linear scaling revealed that the maximum Lyapunov 
exponent survives in the continuum limit, so the real physical system has a 
finite entropy density generation rate. Abelian systems on the other hand 
do not scale linearly leading to vanishing entropy generation in the con- 
tinuum limit. Also analytic studies of simplified subsystems, like the xy 
- model helped to understand the trajectory - defocusing property of the 
potential energy, i.e. the magnetic part of the gauge field interaction 

In order to understand the intriguing coincidence between the numer- 
ically found leading Lyapunov exponent and the gluon damping rate cal- 
culated in high-temperature perturbative gauge theory one has to assume 
that the energy is fully thermalized and the gauge field system evolves to 
a weakly interacting gas of quanta, the gluons. Some arguments in favor 
of a special physical mechanism based on color non-diagonal interaction 
between gluons behind this coincidence has been discussed in [Q, ^ • 

We expect that the non-abelian gauge field system is in its magnetic sec- 
tor the most chaotic. In realistic situations, in particular in ultrarelativistic 
heavy ion collisions there are, however, matter fields: mainly quarks but 
also clusters of several quarks representing a higher color charge. For some 
specific signals, expected to come from such reactions if quark - gluon mat- 
ter has been formed, there are heavy color charges which generate strong 
chromoelectric fields in string - like configurations. As a matter of fact the 
formation and physical behavior of such strings play also a primary role 
in understanding the quark confinement mechanism; the very basis of the 
existence of hadrons. 

It is therefore of genuine interest to investigate the effect of such chro- 
moelectric strings spanned between static charges on the chaotic dynamics 
of random chromomagnetic gauge field configurations. Although specu- 
lative, it also cannot be excluded that we learn about the confinement 
mechanism by simulating the quantum mechanical ground state of QCD, 
i.e. the hadronic vacuum, with random classical chromomagnetic fields. 
Indeed a pair of in total neutral color charges, quantized or classical, con- 
nected by a single fiux line can only be represented by a so called Wilson 
line, the object on which static quark confinement has ever been studied in 
four dimensional euclidean lattice gauge theory. A transition in the chaotic 
dynamical behavior of the classical gauge field system depending on the 
strength of this flux line, on the string constant, may therefore also be 
related to the confinement - deconfinement phase transition. 



2 



Even more speculative but still not totally negligible the idea that ran- 
dom classical field configurations numerically simulate quantum dynamics. 
In fact interference patterns similar to those observed in basic quantum 
mechanical two slit experiments can be created numerically by random 
discrete maps used in certain chaos games 0. Whether the simulation 
reported in the present article has anything to do with such phenomena 
must be a subject of future considerations. 

We organize this paper as follows. First the basic formuli for the clas- 
sical hamiltonian SU(2) lattice gauge field theory are presented. Then the 
speed of energy sharing between chromomagnetic and chromoelectric de- 
grees of freedom with and without static charges is discussed. Finally the 
effect of strong fiux line initialization on the leading Lyapunov exponent 
of the chaotic gauge field system is described. 



In the hamiltonian lattice formulation of classical SU(2) gauge field theory, 
which we use to study its chaotic dynamics, the basic variables are group 
valued on each link of a A^^ cubic lattice]^, 



In the SU(2) model, the smallest nonabelian gauge group, the r^-s are 
the Pauli matrices and the Af{x)-s are corresponding components of the 
nonabelian vector potential. The subscript (x, i) identifies a general lattice 
position X in the interval (0, A^^ — 1) and i the spatial direction of a link 
i = 0,1,2. Since we use periodic lattices this link variable representation 
is complete; the periodic boundary conditions are taken into account in an 
index table of neighboring sites in each of the 6 possible spatial directions. 
The parameters g and a are the bare coupling constant and the lattice 
spacing, respectively. In the classical time evolution actually all space- or 
timelike quantities can be measured in units of a and all types of physical 
energies in units of g^a. 

While the vector potential is related to the basic link variable Ux,i the 
magnetic field strength can be reconstructed from a product of U matrices 
along an elementary closed line, the plaquette: 



In the continuum limit a — > the subleading term is proportional to the 



2 The lattice SU(2) model 
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chromomagnetic field 



U^,ij = exp \^--T''g'^aeijkBlf,j . (3) 

Here €ijk is the totally antisymmetric three dimensional tensor. In practical 
calculations it is advantageous to use the complement link variable, Vx^i, 
instead of the plaquette matrix Up. It is defined so, that the sum of the 
four plaquettes leaning on a given link gives exactly U -V^. 

The chromoelectric field is related to the time derivative of the basic 
link variable, U, 

^". = -ptr(r'^f7.,C/y . (4) 
An exactly good albeit not equivalent definition is 

Again, in practical simulations it is more advantageous to use the canonical 
momentum P = U attached to each link instead. The Hamiltonian of the 
SU(2) lattice gauge theory is given by 

^ = EQ(^'^) + (6) 

in terms of link variables. The summation runs over all links. Here the 
scalar product of two group elements is given by 

{A,B) = ^tv{AB^) (7) 

in the matrix representation and by the usual scalar product of two real 
four dimensional vectors in the quaternion representation. The dynamics 
generated by this Hamiltonian conserves the unitarity of the group ele- 
ments detU =< U,U >= 1, the orthogonality of the basic variables and 
the corresponding canonical conjugates {P, U) = , and commutes with 
the classical Gauss' law. The latter describes the conservation of color 
charge locally, i.e. conserves the following quantity calculated on each 
lattice sites: 

r = ^pc/t-^c/tp, (8) 

+ 

where the "+" summation runs over links originating in the site and the 
"-" over those links which end at the site. 

For the numerical solution of the equations of motion we applied a 
mixed explicit - implicit algorithm which exactly fulfills the constraint of 
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the Noether charge conservation also for discrete timesteps. A derivation 
of this property is given in [10|. The imphcit algorithm, 



U' = U + {P'-eU) 

P' = p + {y-nU + eP') (9) 

with 

e = ([/, P') and = (y, U) + (P', P') , (10) 

where the primed quantities denote the updated (new) values and the el- 
ementary timestep is included in the usage of P = dtU, can be resolved 
into explicit equations. Eventually it can be comprised into the following 
algorithm 

Vt = V - {U,V) U P = P + Vt 

c = ^1 + {P,P) 

P' = {P + U)/c -cU U' = cU + P'. (11) 

It is important to recalculate implicit recursion formuli analytically, be- 
cause explicit equations are not only faster on computer but they also 
conserve the Noether charge (in our case Gauss' law) exactly, i.e. up to 
the highest precision used in the computer code independently of the dis- 
crete time step. An implicit algorithm, like (^, would have been solved 
iteratively so the exactness of the charge conservation depended on the 
iteration depth. 

The explicit numerical algorithm we use conserves Gauss' law, i.e. the 
lattice covariant divergence of the nonabelian electric field (eq.^). The 
only remaining task is to fulfill it initially. For a purely chromomagnetic 
configuration, investigated dynamically up to now, it is trivial to set the 
color charge T on each lattice site to zero by initializing P = on all links. 
This state can also be viewed as a " chromoelectric vacuum". 

To describe localized, i.e. on the lattice pointlike, color charges so that 
the total configuration is color neutral, is far from trivial due to the non- 
abelian nature of Gauss' law in the SU(2) gauge theory. In fact no general 
nonperturbative solution of the classical Gauss' law is known for arbitrary 
color charge density. We follow here, however, a simpler formulation of 
the problem: how to initialize the variables P on the links once a,U - con- 
figuration is given, the total system is color neutral and the charges are 
localized. Our construction principle is as follows. 

1. We generate a C/ - configuration randomly (microcanonically) which 
represents purely chromomagnetic energy. 
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2. We prescribe a flux line by a sequence of 3N neighboring lattice sites. 
In particular we choose a diagonal staircase string across the cube. 

3. On the first link of this flux line Pi = QU\ is initialized. 

4. In order not to have any charge in between we initialize 

Pn = ul_iPn-iUn for (1 < n < 3N) recursively along the prescribed 
flux line. 

5. Arriving at the end of the line the last site has a single incoming 
link with nonzero P on it. Its uncompensated charge is F = —F'^QF 
where F is the path ordered product of basic U variables along the 
flux line (a Wilson line) : 

F = XlUn. (12) 

n 

6. Finally in order to ensure color neutrality of the system we need 
Q — F^QF = for the total charge and tr Q = for ensuring orthog- 
onality of P and U everywhere. The only solution to this problem 
is 

Q = |(Ft-F), (13) 

where the real parameter q describes the strength of the flux line or 
equivalently the magnitude of the end charges. 



3 Less chaos due to flux lines 



We measured the leading Lyapunov exponent as the logarithmic time- 
derivative of the gauge invariant distance, D{t), between two initially ad- 
jacent field configurations 

h=^^\r.D{t). (14) 

This coincides with the original definition of the Lyapunov exponent when- 
ever lnZ)(t) linearly grows. The numerical value of h was fitted in such 
time-intervals only using a time step a = 0.005 equal to the lattice spac- 
ing for the sake of simplicity. The gauge invariant distance of two field 
configurations is defined by 

D{t) = Y.\{U{t),V{t)) - {U*{t),V*{t))\, (15) 

where the superscript * denotes a configuration only slightly different from 
the unstared one initially and the summation runs over all links. This 
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definition of gauge invariant distance differs from the one used earlier |jl] - 
defined using the plaquette sum - only in the actual value of the distance 
but not in its growth rate from which we obtain the leading Lyapunov 
exponent. 

The initially small distance between two trial configurations we produce 
by rotating randomly all link quaternions, U, by an SU(2) group element 
near to unity. Although - in principle - we would have to average the 
value of the leading exponent obtained this way over the phase space, the 
ergodizing property of the chaotic evolution ensures that both trajectories 
sample almost the whole available phase space in a single long term run. 

Fig.|l] shows the evolution of the logarithmic distance. In D{t), for an 
N = 6 cubic lattice with the total scaled energy per link of g^aE = 2.00. 
The different curves counted from above correspond to charges q = 12.10, 
0.0, 8.70, 9.50 at the ends of the flux line and to respectively reduced 
random magnetic backgrounds, belonging to an independent random choice 
of SU(2) group elements with the main group angle covering uniformly 45%, 
100%, 35% and 15% of the interval [0,27r]. 

A tendency of faster initial oscillations in the divergence of the two 
parting configurations with increasing strength of the flux line can be ob- 
served. It is followed by linear epoches with a reduced slope, i.e. reduced 
chaotic dynamics. Correspondingly the entropy generation is also reduced. 

The equipartition of the chromoelectric and chromomagnetic energy 
is essentially faster than the entropy generation rate set by the leading 
Lyapunov exponent. Of course the partition of the total lattice energy 
(not counting the energy of the static charge) develops already in the linear 
approximation without any chaotic behavior. Its time scale is determined 
by the highest oscillation frequency, ~ vr/a, of the lattice system and after 
a few oscillation periods the equipartition is completed. 



A systematic study of the leading Lyapunov exponents as a function 
of the charge strength at the ends of the initial flux line reveals an S-like 
structure known from first order phase transitions (cf. Figj2|). The slight 
initial rise of the upper branch may be due to numerical uncertainties; 
the = 6 system is relatively small in linear size, although the allowed 
phase space is 3239 - dimensional. The qualitative feature of a microcanon- 
ical simulation of a phase transition between weak and strong charges is, 
however, inspiring. 
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Figure 1: The evolution of the gauge invariant distance between two initially 
adjacent configurations with different static charges at the ends of a diagonal fiux 
line in a 6^ cubic lattice. The total energy per link per color is g'^aE — 1.00 for 
all curves. 




Figure 2: The leading Lyapunov exponent as function of the static charges 
positioned at the ends of a diagonal flux line. 



In conclusion these simulations have shown that the chaotic dynamics 
of classical SU(2) Yang-Mills fields survives in the presence of weak static 
charges and a thin flux line of chromoelectric energy spanned between 
them. In order to reduce the chaos essentially a strong (in the order of 
magnitude of 10 in scaled units) charge is necessary. This reduction occurs 
then suddenly at a given critical strength and shows a hysteresis - like in 
case of a ferromagnetic phase transition. 

With respect to relativistic heavy ion collisions this result leaves room 
for the hope that chromoelectric strings (color ropes) formed initially share 
their energy with random chromomagnctic fields - which bundle them to- 
gether - and therefore the soft collective chaotic dynamics of these fields is 
not suppressed essentially. This chaotic dynamics then contributes to the 
thermalization of primordial quark matter helping it to evolve towards a 
quark - gluon plasma. 
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